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I. INTRODUCTION 


Gravitational waves are thought to be produced by astrophysical phenomena ranging 
from the coalescence of orbiting binaries to violent events in the early Universe. Their 
detection would herald a new window for the observation of natural phenomena. Great 
ingenuity is being exercised in attempts to detect such waves in the vicinity of the earth using 
either laser interferometry or various resonant mass devices following Weber’s pioneering 
efforts with aluminum cylinders. Due to the masking effects of competing influences and 
the weakness of gravitation compared with the electromagnetic interactions the threshold 
for the detection of expected gravitationally induced signals remains tantalisingly close to 
the limits set by currently technology. In order to achieve the signal to noise ratios needed 
for the unambiguous detection of gravitational waves numerous alternative strategies are 
also under consideration. These include more sophisticated transducer interfaces, advanced 
filtering techniques and the use of dedicated arrays of antennae. Earth based gravitational 
wave detectors require expensive vibration insulation in order to discriminate the required 
signals from the background. This is one reason why the use of antennae in space offer certain 
advantages. It is argued here that the gravitationally induced elastodynamic vibrations of 
slender material structures in space offer other advantages that do not appear to have been 
considered. Multiple structures of such continua possess attractive properties when used as 
coincidence detectors of gravitational disturbances with a dominant spectral content in the 
10” 4 to 1 Hz region. Furthermore this window can be readily extended to lower frequencies 
and higher sensitivities by enlarging the size of the structures. 

Newtonian elastodynamics [|lj is adequate as a first approximation if supplemented by the 
tidal stresses generated by the presence of spacetime curvature that is small in comparison 
with the size of the detector. The latter are estimated from the accelerations responsible for 
spacetime geodesic deviations. Since the constituents of material media owe their elasticity 
to primarily non-gravitational forces their histories are non-geodesic. The geodesic motions 
of particles offer a reference configuration and the geodesic deviation of neighbours in a 
geodesic reference frame provide accelerations that are additionally resisted in a material 
held together by elastic forces. Since in practical situations the re-radiation of gravitational 
waves is totally negligible the computation of the stresses induced by the tidal tensor of 
a background incident gravitational wave offers a viable means of exploring the dynamical 
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response of a material domain to a fluctuating gravitational field. 

Current resonant mode detectors are designed to permit reconstruction of the direction 
and polarisation of gravitational waves that can excite resonances |2j . Clearly such detectors 
are designed to respond to a narrow spectral window of gravitational radiation and are not 
particularly good at determining the temporal profile of incident gravitational pulses. A 
significant advantage of space-based antennae based on slender material structures is that 
they can be designed to respond to transient gravitational pulses, to polarised uni-directional 
gravitational waves or omni-directional unpolarised waves. 

The general mathematical theory of non-linear Newtonian elasticity is well established. 
The general theory of one-dimensional Newtonian Cosserat continua derived as limits of 
three-dimensional continua can be consulted in [|I| . The theory is fundamentally formulated 
in the Lagrangian picture in which material elements are labelled by s. The behaviour 
of a Cosserat rod at time t may be described in terms of the motion R(s, t ) in space of 
the line of centroids of its cross-sections and elastic deformations about that line. Such a 
rod is modelled mathematically by an elastic space-curve with structure. This structure 
defines the relative orientation of neighbouring cross-sections along the rod. Specifying a 
unit vector d 3 (which may be identified with the normal to the cross-section) at each point 
along the rod enables the state of shear to be related to the angle between this vector and 
the tangent to the space-curve. Specifying a second vector di orthogonal to the first vector 
(thereby placing it in the plane of the cross-section) can be used to encode the state of flexure 
and twist along its length. Thus a field of two mutually orthogonal unit vectors along the 
structure provides three continuous dynamical degrees of freedom that, together with the 
continuous three degrees of freedom describing a space-curve relative to some arbitrary origin 
in space, define a simple Cosserat model (see Figure 1). It is significant for this approach 
that the theory includes thermal variables that can be coupled to the dynamical equations 
of motion, compatible with the laws of thermodynamics. The theory is completed with 
equations that relate the deformation strains of the structure to the elastodynamic forces 
and torques. The simplest constitutive model to consider is based on Kirchoff relations with 
shear deformation and viscoelasticity. Such a Cosserat model provides a well defined six 
dimensional quasi-linear hyperbolic system of (integro-)partial differential equations in two 
independent variables. It may be applied to the study of gravitational wave interactions by 
suitably choosing external body forces f to represent the tidal interaction with each element 
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in the medium. 

A typical system might consist of at least two loops orbiting in interplanetary space. Each 
structure would be composed of high Q material, several km in length. Such a structure, 
made up of transportable segments, could be conveyed to an orbiting station and constructed 
in space. The lowest quadrupole excitation of a steel circular structure would be about 1 Hz 
and vary inversely as its (stress-free) length. Actuator and feedback instrumentation could 
be placed around the antennae to “tune” the system to an optimal reference configuration. 
The precise details of the density and elastic moduli needed to enhance the sensitivity of the 
receiver would result from an in-depth analytic analysis of the Cosserat equations for free 
motion. The ability to readily optimise the resonant behaviour for coupled axial, lateral and 
torsional vibrations by design is a major advantage over other mechanical antennae that 
have been proposed. 

By contrast a broad band detector could consist of an open ended structure coiled into 
a spiral. For planar spirals with traction free open ends the spectral density of normal 
transverse and axial linearised modes increases with the density of the spiral winding number. 
They form ideal broad band detectors with directional characteristics. Furthermore by 
coupling such a spiral at its outer end to a light mass by a short length of high-Q fibre (such 
as sapphire) one may tune such an extension to internal resonances, thereby amplifying 
the spiral clastic excitation. Such excitations offer new detection mechanisms based on the 
enhanced motion of the outer structure of the spiral. 

Cosserat Equations 

The dynamical evolution of a rod with mass density, s G [0, L\ i—>• p(s), and cross-sectional 
area, s G [0, L\ i—>• A(s), is governed by Newton’s dynamical laws: 

pA r = n' + f 
d t (pI(w)) = ml + r' x n + l 

applied to a triad of orthonormal vectors: s G [0, L\ > {di(s, t), d 2 (s, t), d 3 (s, f)} over the 
space-curve: s G [0, L\ » r(s, t ) at time t where n' = d s n , r = d t r, f and l denote external 
force and torque densities respectively and s G [0, L\ i— ► pi is a rod moment of inertia 
tensor. In these field equations the contact forces n and contact torques m are related to 
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the vector w and the strains u, v by constitutive relations. These vectors are themselves 
defined in terms of the configuration variables r and d k for k — 1, 2, 3 by the relations: 

r' — v 


d! k = u x dk 
d k = w x d k . 


The latter ensures that the triad remains orthonormal under evolution, 
identifies 
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Y d k>< d k 

k= 1 


The last equation 


with the local angular velocity vector of the director triad. The general model accommodates 
continua whose characteristics (density, cross-sectional area, rotary inertia) vary with s. For 
a system of two coupled continua with different clastic characteristics on 0 < s < s 0 and 
so < s < L respectively one matches the degrees of freedom at s = so according to a junction 
condition describing the coupling. 

To close the above equations of motion constitutive relations appropriate to the 
rod must be specified: n(s,t) = h(u(s, t), v(s, t), u t (s), v t (s),. .. , s), m(s,t) = 
rii(u(s, t ), v(s, t ), u t (s), v t (s ),... , s) where u t (s) etc., denote the history of u(s , t) up to time 
t. These relations specify a reference configuration (say at t — 0) with strains U(s),V(s) 
such that h(U (s), V(s),... , s) and rh(U(s), V(s),... , s) are specified. A reference configu¬ 
ration free of flexure has r s (s, 0) = d 3 (s, 0), i.e. V(s) = d 3 (s, 0). If a standard configuration 
is such that r(s, 0) is a space-curve with Frenet curvature k 0 and torsion r 0 and the standard 
directors are oriented so that di(s, 0) is the unit normal to the space-curve and d 2 (s, 0) the 
associated unit binormal then U(s) = «o( s ) d 2 (s, 0) + r 0 (s) d 3 (s, 0). 

For a rod of density p and cross-sectional area A in a weak plane gravitational wave 
background the simplest model to consider consists of the Newtonian Cosserat equations 
above with a time dependent body force modelled by a gravitational tidal interaction f. 
In addition to time-dependent waves this may include stationary Newtonian gravitational 
fields. These add terms of the form pAg to f where g is the “effective local acceleration due 
to gravity”. Post Newtonian gravitational fields (such as gravitomagnetic and Lens-Thirring 
effects) can be accommodated with a more refined metric background. 
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An important consideration of any modelling of Cosserat continna to low levels of excita¬ 
tion is the estimation of signal to noise ratios induced by anelasticity and thermal interac¬ 
tions. To gain an insight into the former one may attempt to extend the established theory 
of linear anelasticity to a Cosserat structure. For a string with uniform density p, static 
Young’s modulus E and area of cross section A, the free damped motion in one dimension 
is modelled by the equation: 

pAd u x(s,t ) = n'(s,t) 

where the axial strain v(x,t) = d s x(s,t ) = x'(s,t ) and 

n(s, t) = EA(v(s, t) — 1) — EA f cj)(t — t') v(s,t')dt' 

J — OO 

for some viscoelastic model 0 with 0 < s < L. For free motion in the mode: 

x(s, t ) = s + £(i) cos(7 ts/L) 


the amplitude £(t) satisfies 


€{t) +Wof(t) = / (t>(t - t')£(t')dt’ 


with A? = tt- while for a forced harmonic excitation: 


€(t) + ^(t) =u 0 <j)(t - + F 0 exp(-iQf) 


With ^(0) = ^(0) = 0 the Laplace transformed amplitude of forced axial motion is then 


given in terms of the Laplace transform [1C 


0( a ) = / 0(^) e at dt 


of the anelastic modelling function 0(f) as: 


£(» = 


[a — itt) (a 2 — u>Q<j(f)(<T) + lUq) 


To extend this approach in a simple manner to a 3-D Cosserat model of a slender rod 
with uniform static moduli E and G, geometric elements A, K aa = Jn + J 2 2 , J a p, one 
adopts the following constitutive relations for the local director components of the contact 
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force n and torque m in terms of the local strain vectors v and u and anelastic response 
functions 4>e and 0 G : 


n 3 (s,t) = EA(v 3 (s,t) - 1) - EA / <j> E (t -t')v 3 (s,t') dt 1 


rii(s,t) = GAvi(s,t) — GA / 0 G (f — t') vi(s,t') dt' 


ri2(s,t) = GAv2(s,t) — GA / 0 G (f — t') ^(s, t') dt' 


m 3 (s, t) — K aa G u 3 (s, t) — K aa G / 0 G (f - t') u 3 (s,t') dt 1 


m c 


2 ^ 2 rt 

(s,t) = E'£j 0 #u fl (8,t)-E'EJa P / <t>E(t - t')up(s,t') dt' 
0=1 0=1 


for a, (3 — 1, 2. 


II. BEHAVIOUR OF AN ELASTIC STRING IN A NOISY WEAK GRAVITA¬ 
TIONAL WAVE BACKGROUND 

The above outlines a new approach to the modelling of gravitational interactions with 
slender structures in space. Given the topology and material properties of an antenna one 
may analyse its response to such signals in terms of solutions to a deterministic system of 
well defined partial (integro-)differential equations. These equations are in general easier to 
analyse than those describing the clastodynamics of three-dimensional materials. However 
in addition to controlling the temperature dependence of material characteristics, thermal 
interactions with such structures will induce a stochastic element into the signal response. 
It is therefore necessary to seek modifications to the above Cosserat description that can 
accommodate such random interactions. In the presence of internal damping this is non¬ 
trivial problem for a broad band resonant detector. 

In the calculation that follows we shall oversimplify this problem in order to gain some 
order of magnitude estimates of signal-to-(thermal)noise ratios for both narrow and broad 
band resonant detectors made of known high-Q materials. The simplest approach is to 


7 



approximate the Cosserat equations by ignoring flexure and torsional mode excitations and 
explore the purely string-like excitations. Commensurate with these approximations we shall 
assume that the damping can be described in terms of a single Q-factor at about 4 K for 
each resonant mode and that the thermal interactions give rise to a stochastic system driven 
by a spatially homogeneous Gaussian noise term. 

We shall make two further restrictions by assuming that the antenna possesses sufficient 
stiffness to maintain planar motion in tension (e.g. by overall slow rotation about its centre 
of mass) throughout the excitation time and that the environment produces no overall 
translational drift away from a regular orbital motion of its centroid. The plane will be 
chosen orthogonal to the direction of propagation of a gravitational disturbance. 

We therefore begin by examining the linearised modes about an arbitrary planar 1- 
dimensional structure. In this way we can address the response of both a narrow band 
planar loop and a broad band planar spiral together. Given the anelastic characteristics of 
the structure all of these restrictions can be readily relaxed at the cost of increased com¬ 
plexity in the thermo-mechanical analysis. 

In such a framework consider the (inertial) Cartesian basis {i,j,k} with corresponding 
coordinates (x, y, z ) such that the tidal acceleration field at any position r = xi + yj+zk 
and time t, corresponding to a plane gravitational wave travelling in the direction k at the 
speed of light c has the form jcj 

g(r,t) = z/c)x + A 2 {t-z/c)y^^+ z/c)x-Ai{t-z/c)y^^. (1) 

The dimensionless waveforms A\(t — z/c) and A 2 (t — z/c) arise from two independent polar¬ 
isations of the gravitational wave in the transverse-traceless gauge in linearised gravitation. 
This acceleration excites an elastic string with reference length L , cross-section area A and 
mass density p described by a space-curve r(s,t) parameterised by the parameter s G [0, L\ 
according to 


pAr(s,t)-ri(s,t) = pAg(r(s,t),t). 


The contact force (tension) n(s, t ) is given by the constitutive relation: 


n(s, t) 


EA(V(8,t)\-l) 




( 2 ) 

( 3 ) 
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in terms of the Young’s modulus E. For a closed string, the periodic boundary conditions 
are 

r (0, t) = r(L,t ) 
r'(0,t) = r'(L,t). 

For an open string we consider the tension-free boundary conditions given by 

n(0,t) = n(L,t) =0. (6) 

Let t’o(s) be a stress-free reference configuration of the elastic string with a unit stretch, i.e. 
|t*q(s)| = 1. Assuming the Frenet curvature of r 0 (s), given by 

def 

«o (s) = K(a)| (7) 


(4) 

(5) 


is non-vanishing introduce a Frenet frame as the set of orthonormal vectors along the string: 

i 


def 


to(s) 

= r o(s) 

(8) 

ao(s) 

def *o(s) 

«o(s) 

(9) 


def 


bo(s) 

= t 0 {s) x a 0 (s). 

(10) 


The space-curve of the reference configuration 7*0 (s) is defined in the x-y plane with z — 0 
so that 


With 


r 0 (s) = x 0 (s) i + y 0 (s) j. 


( 11 ) 


b 0 = k (12) 

it follows HH that the tangent vector t 0 (s) may be expressed in terms of the Frenet curvature 
of the reference space-curve as 

t 0 (s) = — sin Kq(s') ds' + d^j i + cos {^J Kq(s') ds' + d^j j (13) 
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for an arbitrary constant 6. Let a(s,t) and /3(s,t) be axial and transverse perturbations in 
the x-y plane of the string about r 0 (s) so that 


r(s,t ) = r 0 (s) + a(s,t)t 0 (s) + /3(s,t) a 0 (s). 

Since the interaction between the string and the gravitational wave is assumed 
shall drop terms containing powers of a, /3, A\, A 2 higher than 1 in the following, 
becomes 

g(r(s,t),t) = ^Ai(t) F^s) + A 2 (t) F 2 (s)J t 0 (s) 

+ | A 2 (t) Fi(s) — Ai(t) F 2 (s)| a 0 (s) 
up to first order terms in a and /?, where 

def 1 

F i(s) = - {^o(s ) 2 -l/o(s) 2 } 

def ^ 

F 2(s) = -{x 0 {s)y 0 {s)}'. 


; From 


, (HD and (HD 


— (1 + OL — Kq /3) F {01 Kq + (3 ,s ) 


a 0 


and therefore 


n' = EA ( a' — kq (5)' t 0 + EA k 0 (a' — k 0 j3) a 0 

up to hrst order terms in a and /3. Substituting this into d|) to first order yields 

E 

a(s,t )-{c/(s, t) — k 0 (s) (3(s, t)Y = F a (s,t ) 

P 

E 

P(s,t) ~ —K 0 (s) {a\s,t) - K 0 (s)P(s,t)} = Fg(s, t) 

where 


def 

Fa(s,t ) = Ai(t) F^s) + A 2 (t) F 2 (s) 


def 

Fp(s,t) = A 2 (t) Fi(s) - Ai(t) F 2 (s). 


(14) 

small we 
Thus (H]) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 
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For an open string the tension-free condition ([|) becomes 


a'(0,t) — n o (0) /3(0,t) = a'(L,t) — k 0 (L) /3(L,t) = 0. (24) 


Introducing 


def 

X(s,t) = a(s,t) 



(25) 


and using (^U[) and fl2T| ) gives 

*(M) = F a (s,t)~ (26) 

which implies 

X{s,t) = .T a (g,t) - +Xo(s) + Xi(s)t (27) 

for two arbitrary functions y 0 ( s ) and Xi( s )- Furthermore by introducing 


f(M) = 


P{s,t) -F 0 (s,t) 
Ln 0 (s) 


(28) 


and substituting (pf5|) and (|27j) into ( ]2T| ) we obtain 

- Av 0 (s) 2 ^(s,t)} = y — {^(s,t) +Xo(s) +Xi{s)t} 

P L P 


(29) 


where 


def 

F(s,t) = ^(s,t) - k 0 (s) ^(sA) = A x (t) T^s) + A 2 (t) F 2 (s) (30) 


in terms of 


def | 

Ti(s) = F[(s) + k 0 (s) F 2 (s) = cos(2/C 0 (s) + 8) 


(31) 


def "] 

F 2 (s) = Fi(s) - K 0 (s) -Fi(s) = -5 sin (2/C„(s) + 8) 


(32) 


with 


def 

/Co (s) = 



Ko(sO ds 7 . 


(33) 
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In deriving the above relations, (|13|) , m and © have been used. It follows from (p5|), 
(p7|) and (|28|) that 


a(s,t) = L£'(s,t) + J 7 a (s,t) + Xo(s) + Xi{s)t (34) 

P(s,t) = Ln 0 (s)£(s : t) + Fp(s,t). (35) 


For an open string with tension-free ends we may substitute (|||) and (|35|) into (|24|), ( |25| ) to 
obtain 


|(0 = 0 


(36) 


by using (|29|) . Therefore a(s,t) and /3(s,t) can be obtained by solving (|29|) for a choice of 
functions yo(s) and Xi(s), subject to the boundary condition (|36l) for an open string or 


^(0 ,t) = £(L,t) (37) 

£'(0 ,t) = e(L,t ) (38) 


together with Xo(0) = Xo(£), Xo(°) = Xo ( L ) and Xi(0) = Xi(L), xi(°) = Xi (L) for a closed 
string. 


III. NORMAL MODE ANALYSIS 


We next analyse ( f20| ) and (^lj) for displacement perturbations that remain small compared 
with L at all times. These solutions represent the deterministic dynamic response of a string 
in terms of small oscillatory deviations from a time-independent reference configuration r 0 (s) 
under the influence of weak gravitational waves. Such a(s, t ) and /3(s, t ) can be obtained by 
choosing 


Xo(s) = Xi(s) = 0 (39) 

that enter (fZ9|), (|34|) and (Po|). For an open string we further choose 

f(0 ,t)=£(L,t) = 0 (40) 

as solutions to (|36|) . To proceed we express £(s, t) in terms of normal modes, a complete set 
of basis functions {^ p (s)} indexed by (integer- or multi-integer valued) p satisfying 

W + = 0 ( 41 ) 
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subject to either periodic or tension-free (/0 P (0) = 0, if) P (L ) = 0) boundary conditions with 
the corresponding eigen-values {-|uIn addition the eigen-functions i/j p (s) shall satisfy 
the ortho-normality conditions 

L 

^p{s) 1 p q (s) ds = L5pq. (42) 

In this basis 



f(M) = JZ&WV’p 


summing over the full range of p. 

Substituting (^9|) and (|43|) into (^) and using ([II]) and (|42|) gives 


£p(t)+vl£ P {t) = f p (t ) 


where 


(43) 


(44) 


def 1 E 

fp(t) = y — ( a i P ^i(f) + «2p^2(f)} (45) 

L p 


with the “overlap coefficients” given by 


def r L 

ai p = — Ei(s) ipp(s) ds (46) 

L Jo 


def r L 

a 2p = - J r 2 {s)^ p {s)ds. (47) 

L Jo 


It follows from (|3T]) and (^) that these coefficients do not depend on the choice of the origin 
on the x-y plane. In addition the coefficients 


def 
dp = 



+ °2p 


(48) 


are independent of 9 and invariant under rotation of the x- and p-axes about the z-axis. The 
value of a p provides a measure of the “coupling strength” between the p-th-normal mode 
and the exciting gravitational wave. 
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IV. MECHANICAL ENERGY IN TERMS OF NORMAL MODES 


The mechanical energy of the string satisfying (|2i) is given by 


£ = 


Substituting (|I4|) into (f49[) gives 


r 2 + EA (|r'| — l) 2 1 ds 


£ = i 
2 


pA ( a 2 + f3 2 ) + EA (a! — /3Y 1 ds 


( 49 ) 


(50) 


up to first order terms in a and (3. By using (|34|) and (|35|) the above expression becomes 


£ = 


J {pA ( L 2 i‘' 2 + L 2 + Tl + Ej-2LiE)+EA (L£" - Lk 2 ^ + rf) d s. 


In the absence of the gravitational wave excitation this reduces to 

2 rL 

i i /-\ i 

° / 0 


P AL I ] H2 


I {?' 2 + «o? 2 + |t 2 Hs 


(51) 


where (^9|) has been used. From (|43|) , ([4l|) and (^2]) together with boundary conditions ( |36| ) 
for an open string or (|37|), ( |38|) for a closed string we obtain 


£ — £ p 


(52) 


with modal contributions 


def 1 / 

£ p - + Q 


(53) 


where 


rrip = 


P 


2 A L 3 u 2 


E 


(54) 


may be identified as an effective modal mass. 


V. SENSITIVITY ESTIMATION 


To accommodate the thermal fluctuations and dissipation for a vibrating string with 


small displacements, the modal equations (|44|) should be modified into a set of coupled 
linear stochastic equations of the form 


€ P {t) + ^22tpqUq£ q {t) +uj 2 p i p {t) = f p (t)+w p (t) 


(55) 
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where the {( pq } are (in general non-diagonal) coupling constants related to material visco¬ 
elasticity and the {w p (t)} are independent Wiener processes characterised by Gaussian prob¬ 
ability densities: 


V(w p (t) - w p (t- At)) = 


exp 


[w p (t) - w p (t-At)Y 
2 b p At 


(56) 


sj2T[bp At 

at any time t for some width parameters b p . As discussed above we assume here that ( pq 
for p ^ q, vanish identically (proportional damping [|(|]). Then (|5^) reduces to a system of 
decoupled linear stochastic equations 


£p{t) A 2 C, p u)p i p (t) + Up £p(t) — f p (t) + w p (t) 


(57) 


def 

in terms of the modal damping ratios ( p = ( pp for any p. These are related to quality factors 


Q p and relaxation times t p by 


Cp — 


2 Q, 


2 OJpTp 


(58) 


To maintain thermal equilibrium at a temperature T between fluctuations and dissipation 
depending on ( p and w p (t) the parameters b p are given by 0 

2 k B T 


bp 


TflpTp 


(59) 


involving Boltzmann’s constant ks (= 1-38 x 10~ 23 J/K.) 

Given initial conditions for £ p (£) and £ p (£) at t = 0, a representative solution of (|57|) for 
t > 0 may be conveniently expressed as 


&(«) = £»(«)+&(«) + &>(«) 


(60) 


where 


eopW 


= fp(0) J 1 - $UpQp(t) + (£ p (0) CpLOp + 4(0)) %{t) 


(61) 


def POO 

«/„(*) = / %(t -t')f r (t')df (62) 

J 0 

def 

Zw P (t) = 


<hp(f — t')w p {t') d t' 


(63) 
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in terms of 


e„(t) = H(t) 


%( t ) = j/(j) 




•P u P 




■p u p 


cos ( \j 1 -Q v p t) 

sin(Jl~(p“pt) 


(64) 


(65) 


using the Heaviside function 


def | 1 if t > 0 

#(*) = 

0 otherwise. 


( 66 ) 


It follows 0 that the £^ p (t) given in (|63|) give rise to the time-dependent deviation values: 

N 2\ k B T 


«» f W 2 > = -S{ i-e-‘ /T -r + } 


(67) 




where 


r ± = 1± sin(2^/l - ( 2 u p t) + ^^sm 2 (^l-Qu p t). 

For given gravitational waveforms Ai(t) and A 2 (t) and the corresponding f p (t ) evaluated 
using (f45|), the “signal” displacements {t) and their derivatives at any measurement time 
t — r > 0 may be calculated using ( |52| ) and compared with those due to thermal “noise” 
given by (0) and (|B8D. 

The above noise terms have small variances if r r p , as pointed out in H in analysing 
a bar-type gravitational wave detector. In this case for “liigh-Q” materials with ( p < 1 
the ratio of “effective mechanical energy” due to signals to “effective thermal energy” at 


<h„m 2 ) = -r-{i-e-‘ /T *r_} 

11 ir) 


K 


( 68 ) 


t = t -C r p may be approximated as follows. Based on 


introduce the effective mechanical 


energy associated with the displacement £j p (t) in (32) due to a signal at measurement time 
r given by: 


def 1 1 

£^{r) = im ^ 


+3^/,(r) 2 


(69) 


= m p £ f (t)f{t)dt-2( p u p mp / £ fp (t) 2 dt. 


(70) 
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Let 


aei 

4 r (i) = 


( 71 ) 


and denote accordingly 




(72) 


so that F v (t) = £,\ t ( t) for 0 < t < r. Now 

JP Jp 

pr poo -| poo 

J o ifp (*) fp(t) dt = J o iff (*) fp(t) dt = ^ 3? J M ^(w)* do; (73) 

in terms of the one-sided Fourier transform operator ~ and its inverse given by 


F(u) = / F(t)e~ lut dt, F(t) = 


2t r 


F(w) e luJt d u 


(74) 


for any function F of non-negative t. It follows from (j72j) that 

where 

1 


$pM = —- 2 . q • 

cn; — u z + 2 z u (Up a; 


and therefore 


1 /*°° 

f/ P (*)/p(*) d * = 7T" 


2 ^p 6c/p 


2tt 7.00 (w 2 - a; 2 ) 2 + 4( 2 a) 2 w 


/!» 


den. 


For sufficiently small (L 


2 C^p OJp LU 


7T 


(o; 2 -n; 2 ) 2 + 4C» 2 


— {5(w — oJp) + i5(u> + Wp)} 


(75) 


(76) 


(77) 


(78) 


£'f\r) 


m r 




— - EALujp 


def ^ 

2 


'p W 'l.p \ 


8P ^p 


l in ([71]) becomes negligible. This yields 

61 Ip (^p) “I - 0^2p ^A-2 (^p) 

(79) 

(|63l) is given by 


2 > + t^p (^p(r) 2 ) 

(80) 

-k B T 

T p 

(81) 
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for r -C r p and ( p C 1. In these limits a gravitational signal to thermal noise ratio can be 
estimated to be 


4 S, (r) 

4 N, M 


u p t p EAL 
2 r ksT 


Clip A-i ((Up) T Cl2p A-2 {u t 


(82) 


The response of any antennae to a general pulse of gravitational radiation can be es¬ 
timated once its spectral content is known. Consider then a linearly polarised harmonic 
gravitational wave with the waveforms: 


A\(t) = h sin u g t 

•A-2 (t) = 0, 


(83) 

(84) 


frequency u 9 /2 ti and dimensionless amplitude h. If the measurement time corresponds to n 

def 

(e Z + ) cycles of this signal i.e. r = r n = 2rm/u g -C t p then 


A^ n (u) = 


hu, 


UJ 2 — Ug 


9 ( e -iwr n _ ^ 


(85) 


Hence ,4i (u p ) has a maximum modulus at resonance when u g = u p , yielding 




ihim 


LU n 


( 86 ) 


In this case the signal-to-noise ratio (^2]) becomes 

£p S \r n ) mrh 2 a\ p Q v EAL 


(87) 


£f>(T„) 4 k„T 

for an integer n -C Q p /2n. This formula is applicable to any 1-dimensional Cosserat structure 
with a specified reference configuration satisfying the material characteristics assumed above. 


VI. A NARROW-BAND CIRCULAR LOOP GRAVITATIONAL ANTENNA 

A circular reference configuration of a closed string may be parameterised with 


ZoO) 

g 

= R cos — 

it 

(88) 

yo(s) 

s 

= R sin — 

it 

(89) 

Zo (s) 

= 0. 

(90) 
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By using (0) this has a Frenet curvature 


(91) 


«o(s) = 


r; 


The “shapes” of the associated normal modes satisfying 
boundary conditions are given by 

V’ImK 5 ) = v/2cos^ 

i>[k, 2 ] (s) = \/2 sin ^ 

IX 


and (|4l|) subject to periodic 

(92) 

(93) 


with multiple mode indices (p = [A:, 1] or [k, 2] with positive integer k). The associated 
eigen-angular frequencies are 


1 / E 

u [k, i] = ^[k, 2 ] = — y (A; 2 + 1) —. 

The effective modal masses, obtained by substituting the above into (U), are 

m[M ]=n»|»,a| = 8ur 3 pyl/? (fc 2 + 1). 


(94) 


(95) 


It follows from (|3l|) and (|32]) with the choice 9 = 0 that 


^ , 1 2s 

•^i(s) = -g cos ^T 

^ , . 1 . 2s 

( s ) = — 2 sm )^- 


(96) 

(97) 


Substituting these into (146| ) and (|47|) we obtain the following overlap coefficients: 

1 


a i 

a 2,[fc,j] = 


2\/2 

1 

2\/2 


«2 ,fc J 
^2,fc ^2J- 


(98) 

(99) 


Therefore only two quadrupole modes (corresponding to p = [2,1] and [2, 2]) are excited by 
the gravitational waves. The signal to thermal noise ratios for these two modes follow from 
(|S2|) as 


A s ) ( T ) 
C [2,j]V r ) 


* u [2j] EAR 
8 r ksT 


AJ (w[2 j]) 


( 100 ) 
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where j = 1,2. From ( |37|) 



( 101 ) 


For R = 1.8 x 10 3 m, p = 8 x 10 3 kg/m 3 , E — 2.0 x 10 11 kg/m s 2 , the quadrupole mode 
frequency is c^[ 2 ,i]/27T = 1 Hz. For a structure with a circular cross-section of radius r = 0.01 
m the area A = nr 2 . If one selects T — 4 K, Q[ 2 ,i] = 10 5 and n ps 0.1 x Q[ 2 ,i]/27t then the 


condition £^\] ( r n)/( r n) = 1 hnplies an amplitude sensitivity of h ~ 2 x 10 21 for such 


an antenna. 

VII. A BROAD-BAND SPIRAL GRAVITATIONAL ANTENNA 

By contrast to the narrow-band loop antenna in which only two quadrupole modes with 
the same frequency are excited by a harmonic gravitation wave as discussed in Section |Vl]. 
we now demonstrate how a multi-frequency gravitational antenna may be constructed based 
on a spiral reference configuration. Consider an open string whose reference configuration 
lies on the x-y plane as before but with a non-constant Frenet curvature of the form 



( 102 ) 


for a constant p and scale parameter R. It follows from (|33|) that 


2 P 3 

Ms) = o— IMs ) 3 - «o(o) 3 ]. 

o fJj 


(103) 


Substituting ( |102| ) into ([4d]) yields 



(104) 


where 



(105) 


Solving ( |104|) subject to tension-free boundary conditions yields 



( 106 ) 
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and the condition 


Ai(-Ap) Bi(-Ap) 


A p 


0 

Ai(p£ — Ap) Bi(pf — Ap) 


B p 


0 


( 107 ) 


def 

where £ = L/R and the Ai, Bi are the standard Airy functions. The corresponding charac¬ 
teristic equation for X p , 

Ai(— X p ) Bi(pf — X p ) — Ai (fj,£ — X p ) Bi(— X p ) = 0, (108) 


determines the eigen-values A p . When comparing the spectrum of the eigen-frequencies of a 
spiral with that of a circular loop it is convenient to introduce the function 


fi(s) 


def 


«o(s) 



(109) 


Note that for any s G [0, L] the spiral has a radius of curvature 1/ko(s) and from (|94|) it 
can be seen that f2(s) equals the quadrupole mode angular frequency of a circular loop of 
radius 1/k,q(s) with the same mass density p and Young’s modulus E. Given values of X p 
as solutions of (|108| ), it follows from ( |105| ) and ( |109| ) that the eigen-angular frequencies u p 
are given by 


u p = n(o) 


p? X p +1 


( 110 ) 


and from (|107|) 


o Ai(/A' - A p ) A 
p _ Bi(p£ - A p ) 

The normalisation condition (0)3) now takes the form 


( 111 ) 


{A p Ai (pa — A p ) + Bi (per — A p )} 2 dcr = £ 


( 112 ) 


where a = s/R. Therefore ( |111| ) and ( |112| ) determine A p and B p . Furthermore the overlap 
coefficients ai p and d 2 P are obtained from (^6j) and (|47| ) using ([JT]), (^2|) and 
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Finally from flS7| ) the signal to (thermal) noise ratio for the p-th mode of the spiral may 
be written: 

4 5) 0n) nir£h 2 a 2 lp Q p EAR 

Si N \r n ) ~ ^ k B T [ 

In Figure 4 we have plotted, for a given spiral geometry (determined by the parameters 
/i, i and R, the values of a 2 = ^(a^ p +a| p ) for the first 60 normal modes of the spiral according 
to the equations above. The modes are distributed uniformly on the abscissa in terms of the 
non-dimensional eigen-angular frequencies cD p = u p /Q(0). The broad-band response feature 
of this antenna is clearly visible in the region where significant signal overlap occurs. 

For /i = 0.3, £ — 167r, R — 1.9 x 10 3 m, p = 8 x 10 3 kg/m 3 , E = 2.0 x 10 11 kg/rn s 2 , 
one finds for p = 45 that X p = 95.5 and a lp = 0.135 corresponding to the eigen-frequency 
u p /2i t = 1 Hz. For a structure with a circular cross-section of radius r = 0.01 m the area 
A = 7ir 2 . If one selects T = 4 K, Q p = 10 5 and n ~ 0.1 x Q p /2tt then the condition 
4 S) (r n )/(r n ) = 1 implies an amplitude sensitivity of /r ~ 2 x 1CT 21 for such an antenna. 
Thus the spiral antenna maintains a sensitivity commensurate with that of the single loop 
antenna but with the added broad-band characteristic. 

VIII. CONCLUSIONS 

The simple Cosserat theory of rods has been outlined and used to estimate the thermal 
noise sensitivity of gravitational antennae constructed out of orbiting slender material struc¬ 
tures. Although such calculations have not included many other competing noise sources 
we feel that they provide support for a novel concept. Orbiting Cosserat structures can 
accommodate both narrow-band and broad-band detectors and may offer a much cheaper 
alternative to existing space-based laser-interferometers. One may venture optimism that 
costs will become ever more competitive with the current rapid development of high-strength 
carbon-based fibres. Not withstanding economic considerations, the implementation of these 
ideas would be complimentary to existing global programmes of gravitational wave research. 
The design of detectors with optimised response in the 1 Hz spectral region would exploit a 
current niche in the existing gravitational antenna frequency spectrum and the detection of 
gravitational waves in this domain would provide vital information about stochastic back¬ 
grounds in the early Universe and the relevance of super-massive black holes to processes in 
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the centre of galaxies. 

Although we have concentrated on planar structures and short measurement times, in 
order to gain enhanced sensitivities, the extension to non-planar loops and spirals is in 
principle straightforward. For example, several “wire-balls” of approximately spherical shape 
may offer a viable method to monitor stochastic gravitational waves over longer periods of 
time. 

There remain many further important issues that need detailed attention. However we 
believe that the use of stochastic Cosserat thermo-mechanics for slender rods in noisy gravi¬ 
tational fields offers a reliable research tool for the development of the scenarios outlined in 
this paper. 
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FIG. 2: Image of a spiral with fi = 0.3 and t = 16 7r plotted in the - y) plane. The radius of 
curvature of this spiral varies from l/«;o(0) = R at the outer end to \/kq{L) = 0.651? at the inner 
end with approximately 8 (~ 1/2-k) windings. 
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FIG. 3: Behaviour of the first 3 eigen-functions, V’l(s) (solid), ip 2 (s) (dotted) and -03 (s) (dashed) 
for a spiral with fi = 0.3 and i = 16 7r, plotted against o = s/R 6 [0,£]). 
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FIG. 4: The overlap parameters a p = £{a\ p + a^) for the first 60 normal modes (1 < p < 60) 
of a spiral with /i = 0.3 and i = 16 7r plotted against the non-dimensional angular frequency 
Cb p = u; p /O(0) calculated from \ p using (11T). Relatively larger a p are obtained for 30 < p < 50 
corresponding to 1 < tu p < kq(£)=1.54: between the two vertical lines. 
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